© BUNDESREPUBLIK 
DE UTS CH LAND 




DEUTSCHES 
PATENTAMT 



© Pat ntschrift 
© DE 19505 283 C1 



© Aktenzeichen: 195 05 283.6-33 

© Anmeldetag: 16. 2. 95 

© Offeniegungstag: — 
© Veroffentlichungstag 

der Patenterteilung: 11. 4.98 



g) IntCI. 8 : 

H05G 1/64 

G 01 N 23/04 
H 04 N 5/32 
A81B6/00 



S 

CM 
1A 

U> 



Innerhalb von 3 Monaten nach Veroffentlichung der Erteilung kann Einspruch erhoben warden 



LU 

a 



© Patentinhaber: 


© Erfinder: 


Siemens AG, 80333 Munchen, DE 


Hassler, Dietrich, Dipl.-lng., 91080 Uttenreuth, DE 




© Fur die Beurteilung der Patentfahigkeit 




in Betracht gazogene Druckschriften: 




EP 04 87 108 A1 




»NucL instr. and Methods in Physics Research* A 




310(1991) 460-484; 



© Rontgenuntersuchungsanlage 

© Es soil sine Rontgenuntersuchungsanlage mit zwei bildge- 
benden Systemen geschaffen warden, bei der Bildartefakte 
durch Straustrahlung weitgehend reduziert sind. 
Die bildgebenden Systeme welsen a-Si:H-Detektoren auf, 
welche wechselsweise mit Rontgenpulsen (4) beaufschfagt 
warden, zwischen denen die Bildauslesung erfolgt, wobei 
der R5ntgenpul8 (4) des einen Detektors in der Pause (7) 
zwischen der Auslesung der beiden Halbbilder des anderen 
Detektors liegt und der Storanteil durch Drfferenzbildung 
zweier Halbbilder ermlttelt und vom Gesamtbild abgezogen 
wird. 
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In der Rontgenkardiologie werden zwei gleichzeitig 
aktive, gekreuzte, bildgebende Systeme (Zweiebenen- 
anlage) betrieben. Die unvermeidliche Streustrahlung 
aus dem Patienten beirn Betrieb des einen Systems" er- 
reicht auch den Detektor des anderen Systems und f Qhrt 
dort zu unerwiinschten Oberlagerungen im Bild. Dies 
gilt wechselseitig. Eine Zweiebenenanlage ist beispiels- 
weise in EP 0 487 108 Al beschrieben. Dab ei sind CCD- 
Detektoren vorgesehen. Damit die Streustrahlung nicht 
zu einer Verfaischung des Bildes fuhrt, wird der CCD- 
Detektor vor dem Auftreten des Streustrahlungsimpul- 
ses in einen far die Bildaufnahme inaktiven Zustand 
versetzt In "NucL Instr. and Methods in Physics Rese- 
arch" A310 (1991), S. 460—464 ist die Verwendung von 
a-Si:H-Detektoren in bildgebenden Systemen beschrie- 
ben. Bei dieser Art von Detektoren ist die Unterdriik- 
kung von Streustrahlungsartefakten in der geschilder- 
ten Weise nicht mSglich. 

Der Erfindung liegt die Aufgabe zugrunde, eine R6nt- 
genuntersuchungsanlage mit zwei bildgebenden Syste- 
men mit a-Si:H-Detektoren so auszubilden, daB die 
durch Streustrahlung erzeugten Artefakte weitgehend 
reduziert sind. 

Diese Aufgabe ist erfindungsgemaB geldst durch die 
Merkmale des Patent anspruchs 1 . 

Weiterbildungen und Einzelheiten der Erfindung er- 
geben sich aus dem Unteranspruch. 

Die Erfindung ist nachfolgend anhand der Zeichnun- 
gen naher erlautert Es zeigen: 

Fig. 1 Pulsverlaufe zur Erlauterung des Erfindungsge- 
dankens, 

Fig, 2 eine schematische Darstellung des Detektors 
eines bildgebenden Systems fur eine Rontgenuntersu- 
chungsanlage nach der Erfindung, 

Fig, 3 und 4 zwei Varianten zur Erlauterung der erfin- 
dungsgem&Ben Artefaktreduktion, und 

Fig. 5 bis 18 Kurvenverlaufe zur erfindungsgemaBen 
Artefaktreduktion. 

In der Fig. 1 sind in der oberen Zeile die Pulse eines 
ersten bildgebenden Systems dargestellt Ein Puis 1 ver- 
einigt einen Teilpuls 2 zur Loschung der Restinforma- 
tion vora vorhergehenden Bild, einen Teilpuls 3 zur Auf- 
ladung der Photodioden in dem a-Si:H-Detektor des 
bildgebenden Systems und einen Rdntgenpuls 4. Zwi- 
schen zwei Pulsen 1 liegt die Auslesezeit, die fur das 
erste Halbbild mit 5 und das zweite Halbbild mit 6 be- 
zeichnet ist. Dazwischen liegt eine Pause 7. 

In der unteren Zeile der Fig. 1 wird gezeigt, daB die 
Pulse 1 des zweiten bildgebenden Systems zeitlich ver- 
setzt gegenuber den Pulsen 1 des ersten bildgebenden 
Systems sind. Ein Rdntgenpuls 4 des zweiten bildgeben- 
den Systems erzeugt eine unerwiinschte Streustrahlung, 
die auf den Detektor des ersten Systems falit und dort 
auf die noch nicht ausgelesene Biidinformation addiert 
wird. Bei ununterbrochenem Betrieb ist diese Stoning 
nicht vermeidbar, kann aber kompensiert werden, wenn 
folgendes gilt: 

a) Die Streustrahlung allein fuhrt zu einem Bild, das 
im Vergleich zur Systemgrenze, gegeben durch den 
Zeilenabstand, nur tiefe vertikale Ortsfrequenzen 
(max. 1/2 Nyquist-Netzfrequenzband = 1/4 Ab- 
tastfrequenz) enthalt 

b) Die Systeme werden so synchronisiert, daB die 
Pulse 1 jeweils nach der Halfte der Auslesezeit des 
anderen Systems kommen (Fig. 1). 



c) Das Bild kann im Zwischenzeilenverfahren (In- 
terlaced Mode) ausgelesen werden, d. h. zum Bci- 
spiel zuerst in der ersten Halfte der Auslesezeit ein 
Halbbild aus nur den ungeraden Bildzeilen, dann 
nurdengeraden. 

d) Wahrend der Strahlungszeit des jeweils anderen 
Systems, zwischen der Auslesung von zwei Halbbil- 
dern, wird eine Pause der Auslesung eingeschaltet 
(Fig-1). 



10 



Das erste Halbbild ist von storenden Oberlagerungen 
frei, weil der Puis 1 des anderen Systems noch nicht 
auftrat. Das zweite Halbbild enthalt das Nutzbild mit 
Oberlagerter Storung. Die Differenz beider Bilder ent- 
15 halt fur tiefe Frequenzen (max. Frequenz < 1/4 y-Ab- 
tastfrequenz) den Storanteil allein, Er kann vom Ge- 
samtbild abgezogen werden. 

Die Fig. 2 zeigt schematisch den a-Si:H-Detektor 5 
eines bildgebenden Systems mit den ungeraden Zeilen 6 
20 und den geraden Zeilen 7. 

Den KorrekturprozeB zeigt Fig. 3. In Fig. 3 werden 
auf der Leitung 8 die den geraden Zeilen 7 entsprechen- 
den Bildsignale und auf der Leitung 9 die den ungeraden 
Zeilen 6 entsprechenden Bildsignale TiefpaBfilter 10, 11 
25 zugefuhrt deren Ausgangssignale nach Interpolation 
auf die voile Zeilenzahl in den Interpolatoren 12, 13 
einer Subtraktionsstufe 14 zugefuhrt werden. Deren 
Ausgangssignal wird zusammen mit den Signalen auf 
den Leitungen 8, 9 einer Additionsstufe 15 zugefuhrt, 
30 deren Ausgangssignal am Ausgang 16 ein von Streu- 
strahienartefakten weitgehend freies Bild verkorpert. 

Aus den beiden Halbbildern werden durch HefpaBf il- 
terung (Sine-Interpolation) in y-Richtung (zweidimen- 
sionale Orts-TiefpaBfilterung fuhrt zu Ergebnissen, die 
35 von Ungleichheiten der Kanale unabhangiger sind) und 
Interpolation der fehlenden Zwischenzeilen Bilder der 
vollen Zeilenzahl erzeugt Nach Subtraktion liegt das 
Streubild isoliert vor. Die Nutzinformation, die in bei- 
den Teilbildern gleich ist, wird zu Null subtrahiert, wenn 
40 Unsymmetrien der Kanale vermieden werden. Das 
Streubild wird dann vom gestdrten Vollbild abgezogen. 
Lediglich hochfrequente Anteile des Quantenrauschens 
des Streubildes (jenseits der TiefpaB-Eckf requenz) blei- 
ben als Reststorung unkompensiert zuriick. 
45 Filterung und Interpolation werden vorzugsweise di- 
gital als Faltung mit einem Faltungskern durchgefuhrt, 
der Stutzwerte fur die voile Zeilenzahl enthalt 

Eine vereinfachte Version nach Fig. 4 ist mdglich, 
wenn man Fehler vorzugsweise im Frequenzbereich na- 
50 he der Nyquistgrenze der Halbbilder zulaBt Nach zum 
Beispiel linearer Interpolation aus Nachbarwerten ge- 
nugt dann eine einzige TiefpaBfilterung im Filter 17 im 
Differenzkanal am Vollbild. Liegt die Eckfrequenz des 
Filters 17 niedrig, d. h. < 1/4 y-Abtastfrequenz, werden 
55 die Fehler unterdrtickt 

Tiefe Ortsfrequenzen in y-Richtung mQssen in beiden 
Halbbildern trotz des Zeilenversatzes praktisch gleich 
wiedergegeben werden. DaB dies angenommen werden 
darf, zeigt die folgende Systembetrachtung: 
60 Die Fig. 5 und 6 behandeln zunachst die Signalverarbei- 
tung im normalen Bild voller Zeilenzahl. Dem Ortsraum 
in Fig. 5 ist immer der Fourier-Raum in Fig. 6 zugeord- 
net Eine sinusformige Helligkeitsverteilung in y-Rich- 
tung, quer zu den Zeilen, ist Ausgangspunkt in Zeile a. 
65 Sie wird raumlich beschnitten (endliche DetektorgroBe; 
Zeile b), was zum Spektrum in Zeile c (Fig. 6) fuhrt Die 
Abtastung an acht Stellen pro Periode im Abstand sp 
(Zeile d; Abtastf requenz 1/sp) fuhrt auf das Spektrum in 
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Zeile e. Die Abtastung filhrt zu hdherfrequenten Spek- 
tralanteilen, die zur eindeutigen Rekonstruktion des ur- 
sprfinglichen Signalverlaufs (Zeile g) mit einem TiefpaB 
(Zeile f) hdchsteris der Eckfrequenz 1/2 sp (Nyquist-Fre- 
quenz) weggefiltert werden k6nnen. Der Filterung (un- 5 
endlicher Steilheit) entspricht eine Faltung mit der sinc- 
Funktion (Faltungskern). 

Dieselbe Helligkeitsschwankung niedriger Frequenz 
wird nun in den Fig. 7 und 8 mit halbierter Zeilenzahl, 
z. B. mit alien ungeraden Zeilen, abgetastet Das Prinzip 10 
ist unver&ndert. Nur die Abtast- und Nyquist-Frequenz 
sind halbiert Die Rekonstruktion des Originals aus der 
abgetasteten Funktion durch Tiefpasse ist noch genauso 
m&glich, weil die Spektren sich im Beispiel noch nicht 
uberlappen und die Grenzfrequenz des Tiefpasses eben- 15 
falls halbiert ist, so daB die Trennung der Spektralanteile 
gelingt Im Grenzfall unendlich ausgedehnter Eingangs- 
schwingungen ist das Spektrum der Eingangsschwin- 
gung unendlich schmal und die rekonstruierbare Fre- 
quenz kann bis zur halbierten Nyquist-Grenze reichen. 20 
Eine praktische Grenze fur die TiefpaB- Eckfrequenz 
liegt bei ca.3/4 dieses theoretischen Grenzwertes. 

Die Fig. 9 und 10 zeigen die Abtastung nur mit den 
geradzahligen Zeilen. Obwohl die Abtastung zu ande- 
ren Werten fuhrt (Zeile e), ergibt die Rekonstruktion 25 
der Eingangsschwingung durch TiefpaBfilterung/Fal- 
tung dasselbe Ergebnis wie in den Fig. 8 und 9. 

Damit ist belegt, daS die tieffrequenten Spektraliantei- 
le des Nutzbildes — bis auf Ungleichheiten in den Kana- 
len — durch die Subtraktion nach Interpolation (Fig. 3) 30 
herausfallen und nur der niederfrequente Streuanteil 
unverfaischt erhalten bleibt, so daB er vom Gesamtbild 
abgezogen werden kann. 

Die hdherfrequenten Spektralanteile (Nutzbild und 
Quantenrauschen des Streubildes) fuhren aber zu un- 35 
vermeidlichen Artefakten. Da der vor dem Detektor 5 
liegende Szintillator in seiner Eigenschaft als Anti-Alia- 
sing-Filter des Eingangsbildes auf die voile Zeilenzahl 
abgestimmt ist, kommen Frequenzen zwischen 1/4 und 
1/2 der Abtastfrequenz (1-Pixelabstand) vor und fuhren 40 
zu Doppeldeutigkeiten in den Teilbildem. Dies ist in den 
Fig. 1 1 bis 1 4 dargestellt 

Die Frequenz des Eingangssignals ist im Beispiel 3/4 
der Abtastfrequenz der Teilbilder. Dies fuhrt in den 
Fig. 1 1 und 12 trotzdem zu gleichen Abtastwerten (Zeile 45 
e) wie in den Fig. 7 und 8. In den Fig. 13 und 14 sind die 
Abtastwerte invertiert. Im Spektrum (Zeile e) werden 
die hochfrequenten Eingangssignale in tiefe Lagen hin- 
abgemischt (Aliasing). Die Phasenumkehr (Zeile g) fuhrt 
dazu, daB sich die Beitrage in der Subtraktion nicht 50 
auslSschen, sondern sogar addieren. 

Eine Minderung dieses unerwunschten Nebeneffek- 
tes kann man auf folgende Weise anstreben. Pro Voll- 
bild wird zwischen ungerad- und geradzeiligen Teilbil- 
dem fur die Auslesezeiten vor/nach dem Storpuls des 55 
fremden Systems gewechselt Damit wechselt die Pola- 
ritat der Aliasing-Storung von Bild zu Bild. Wegen der 
Tragheit des Auges wird dies zu einer Kompensation 
oder Verwischung der Storung fuhren, wenn die Bildfre- 
quenz hoch genug und die Bewegung des Objektes nicht 60 
zu schnell ist Staranteile aufgrund von Rauschen sind 
nicht kompensierbar. Sie addieren sich quadratisch un- 
ter der Wurzel. 

Da der angesprochene Frequenzbereich bei typi- 
schen MTFs des Szintillators nur schwach vertreten ist, 65 
wird davon ausgegangen, daB die Aliasing-Artefakte — 
bereits ohne den obengenannten Polaritatswechsel — 
weniger sttirend auffailen als die Streustrahlung vom 



Nachbarsystenx 

Zur Erklarung der vereihfachten Signalverarbeitung 
nach Fig. 4 dienen die Fig. 15 bis 18. Hier ist eine noch 
niedrigere Eingangsfrequenz und eiri lingerer Aus- 
schnitt angenommen. Die zu trennenden Spektralkom- 
ponenten (Zeile e) liegen so weit auseinander, daB sogar 
ein weniger steiles TiefpaBfiher ausreicht Eine einfache 
Realisierung besteht z. B. in der linearen Interpolation 
von Zwischenwerten aus den vorhandenen Nachbar- 
werten. Man erkennt dies auch an den Abtastwerten in 
Zeile e, die sich bis auf den Phasenunterschied 1 /sp 
kaum unterscheiden. 

Bei steigender Eingangsfrequenz werden die Fehler 
immer groBer. Da aber auch hochfrequentere Nutzsi- 
gnale in beiden Teilbildem gleichartig erfaBt werden 
mussen, um sich bei Subtraktion zu kompensieren, diir- 
fen sie nicht zur Korrektur zugelassen werden. Daher 
muB ein TiefpaB geringer Eckfrequenz gewahlt werden. 

Da die Streustrahlungsbilder sehr niederfrequent 
sind, geniigt eine grobere Abtastung. Statt jeder zwei- 
ten Zeile genugt es, jede 16. Zeile abzutasten, um das 
Streustrahlungsbild noch hinreichend gut zu erfassen. 
Dementsprechend kurzer ist die zugeordnete Auslese- 
zeit. Die Zeilenzahl des Teilbildes ohne Storung ist ent- 
sprechend hoher und seine Auslesezeit verl&ngert Der 
Rontgenpuls des zweiten Systems rilckt damit naher an 
den nachsten Puis des ersten Systems. Die Gleichzeitig- 
keit der Rontgenpulse ist besser angenahert. 

Patentanspriiche 

1. Rontgenuntersuchungsanlage mit zwei bildge- 
benden Systemen mit a-Si:H-Detektoren (5), wel- 
che wechselweise mit Rontgenpulsen (4) beauf- 
schlagt werden, zwischen denen die Bildauslesung 
erfolgt, wobei der Rontgenpuls (4) des einen Detek- 
tors (5) in der Pause (7) zwischen der Auslesung der 
beiden Halbbilder des anderen Detektors (5) liegt 
und der Stdranteil durch Differenzbildung zweier 
Halbbilder ermittelt und von den Zeilen des Ge- 
samtbilds abgezogen wird, die Streustrahlenanteile 
enthalten. 

2. Rontgenuntersuchungsanlage nach.Anspruch 1, 
bei der aus den beiden Halbbildern durch TiefpaB- 
filterung und Interpolation der fehlenden Zwi- 
schenzeilen Bilder der vollen Zeilenzahl erzeugt 
werden, welche voneinander subtrahiert werden. 
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